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We study the relation between the spectral gap above the ground state and the decay of 
the correlations in the ground state in quantum spin and fermion systems with short-range 
interactions on a wide class of lattices. We prove that, if two observables anticommute 
with each other at large distance, then the nonvanishing spectral gap implies exponential 
decay of the corresponding correlation. When two observables commute with each other 
at large distance, the connected correlation function decays exponentially under the gap 
assumption. If the observables behave as a vector under the U(l) rotation of a global sym- 
metry of the system, we use previous results on the large distance decay of the correlation 
function to show the stronger statement that the correlation function itself, rather than 
just the connected correlation function, decays exponentially under the gap assumption on 
a lattice with a certain self-similarity in (fractal) dimensions D < 2. In particular, if the 
system is translationally invariant in one of the spatial directions, then this self-similarity 
condition is automatically satisfied. We also treat systems with long-range, power-law 
decaying interactions. 
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1 Introduction 

In non-rclativistic quantum many-body systems, a folk theorem states that a nonvanishing 
spectral gap above the ground state implies exponentially decaying correlations in the 
ground state. Perhaps this has been the most popular folk theorem in this field since 
Haldane [1] predicted a "massive phase" in low dimensional, isotropic quantum systems. 
Quite recently, this statement was partially proved [2] for quantum lattice systems with a 
global U(l) symmetry in (fractal) dimensions D < 2. More precisely, abound which decays 
to zero at large distance was obtained for correlation functions whose observables behave 
as a vector under the U(l)-rotation. Unfortunately, the bound is weaker than the expected 
exponential decay. On the other hand, exponential clustering of the correlations was also 
proved recently [3, 4] for quantum many-body lattice systems under the gap assumption. 
This is a non-relativistic version of Fredenhagen's theorem [5, 6] of relativistic quantum 
field theory. Clearly the following natural question arises: can this clustering property 
be combined with the above bound for the decay of the correlations to yield the tighter, 
exponentially decaying bound for the correlation functions themselves, rather than just 
for the connected correlation functions? We emphasize that these are different statements; 
given clustering, the decay of the correlation functions requires also that certain matrix 
elements vanish in the ground state sector. 

In this paper, we address this problem and reexamine the above folk theorem by relying 
on the exponential clustering of the correlations. Our first step is to provide a rigorous 
proof of the exponential clustering. We extend the previous results in this case to treat 
long-range interactions including both power-law and exponentially decaying interactions. 
In the former case, all the upper bounds for the correlations become power-law bounds. 

We then prove that ground state correlation functions of observables which transform 
as vectors under a U(l) symmetry decay exponentially or with a power law, depending on 
the form of the interaction, given an additional assumption on a certain self-similarity. In 
particular, if the system is translationally invariant in one of the spatial directions, this 
self-similarity condition is automatically satisfied. Therefore the corresponding correlation 
functions decay exponentially for translationally invariant systems on one-dimensional 
regular lattices. As a byproduct, we also prove that, if two observables anticommute with 
each other at large distance, then the corresponding correlation in the ground state decays 
exponentially imdcr the gap assumption for a wide class of lattice fermion systems with 
exponentially decaying interactions in any dimensions. In this case, we do not need any 
other assumption except for those on the interactions and the spectral gap. 

This paper is organized as follows: In the next section, we give the precise definitions of 
the models, and describe our main results. In Section 3, we prove the clustering of generic 
correlation functions under the gap assumption, and obtain the upper decaying bound for 
the fermionic correlations. The decay of the bosonic correlations are treated in Section 4. 
Appendix A is devoted to the proof of the Lieb-Robinson bound for the group velocity 
of the information propagation in the models with a long-range interaction decaying by 
power law. 
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2 Models and main results 

We consider quantum systems on generic lattices [7]. Let be a set of the sites, 
x,y, z,w, . . ., and a set of the bonds, i.e., pairs of sites, {x,y}, {z,w}, . . .. We call 
the pair, A := (As,Ab), the lattice. If a sequence of sites, Xo,Xi,X2,.--,Xn, satisfies 
{xj-i,Xj} e A(, for j — l,2,...,n, then we say that the path, {xo,xi,X2, ■ ■ ■ ,Xn}, has 
length n and connects xq to Xn- We denote by dist{x,y) the graph-theoretic distance 
which is defined to be the shortest path length that one needs to connect x to y. We 
denote by |X| the cardinality of the finite set X. The Hamiltonian H\ is defined on the 
tensor product <S>x&As of a finite dimensional Hilbert space Tlx at each site x. We 
assume sup^^ sup^. dim Tlx < N < oo. For a lattice fermion system, we consider the Fock 
space. 

Consider the Hamiltonian of the form, 

Ha^ Yl hx. (2-1) 

XcAo 

where hx is the local Hamiltonian of the compact support X. We consider both power-law 

and exponentially decaying interactions hx- 

For the power-law decaying interactions hx, we require the following conditions: 

Assumption 2.1 The interaction hx satisfies 

with positive constants, Aq and rj, and the lattice A equipped with the metric satisfies 

V - X - < ^ (2 3) 

f^^^ [1 + dist(a:, z)\^ [1 + dist(^, y)Y " [1 + dist(x, y)\^ ^ " ' 

with a positive constant po- 

Remcirk: If ^ 

sup sup — -. — < oo, (2.4) 

then the inequality (2.3) holds as follows: 

}2 ^- — ^- — ^ X 



< 



^^^^ [1 + dist(x, ^)]'' [l + dist{z,y)]^ 

[1 + dist(x, y)]" [1 + dist(,T, ^)]''[1 + dist(^, y)]i 

1 ^ [1 + dist(x, z)Y' + [1 + dist(^, y)]'' 

[1 + dist(a;, y)]"^ [1 + dist(a;, z)]^[l + dist(;s, y)]"^ 



^ ]_ 2^ I ^ I ^ I (2 5) 

- [l + dist{x,y)]\^^ \ [1 + dist(x, ^)]'? [l + dist(^,7/)]^j ' ^'^ 
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where we have used the inequahty, [1 + dist(a;, y)]'^ < 2^([1 + dist(a;, z)]'^ + [1 + dist(z, y)]^). 
Prom the assumption (2.2) and the condition (2.4), one has 

sup ^ < So < oo, (2.6) 

where sq is a positive constant which is independent of the volume of |As|. 
Instead of these conditions, we can also require: 

Assumption 2.2 The interaction hx satisfies 

sup J2 ||/ix|||X|[l + diam(X)]'' < Si < oo, (2.7) 

^ X3X 

where rj is a positive constant, diam(X) is the diameter of the set X, i.e., diam(X) = 
max{dist(x, y)| X, y e X}, and Si is a positive constant which is independent of the volume 
of\As\- 

For exponentially decaying interactions hx, we require one of the following two as- 
sumptions: 

Assumption 2.3 There exists a positive rj satisfying the condition (2.4)- The interaction 
hx satisfies 

E \\hx\\<\oexp[-{fi + e)dist{x,y)] (2.8) 

X3x,y 

with some positive constants, Xq, /i ande. 
Remark: From the conditions, we have 

r / N -, / M An exp[— /i distfx. v)l 
exp|-(^ + .)d.t(x.,)]< (2.9) 

with a positive constant Aq, and 

^ exp[-/idist(x,z)] ^ exp[-/x dist{z,y)] ^ poexp[-/i dist{x,y)] 

[l + dist(a;,2;)]'? ^ [1 + dist(2;, 2/)]" " [1 + dist(a;, 2/)]'? ^ ' ' 

with a positive constant po in the same way as in the preceding remark. 



Assumption 2.4 The interaction hx satisfies 

sup E exp[/i diam(X)] < si < oo, (2-11) 

^ X3X 

where is a positive constant, and Si is a positive constant which is independent of the 
volume of \As\. 
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Remcirk: This assumption is milder than that in [6] by the absence of the factor A^^l^l 
in the summand. 

Further we assume the existence of a "uniform gap" above the ground state sector of 
the Hamiltonian H\. The precise definition of the "uniform gap" is: 

Definition 2.5 (Uniform gap): We say that there is a uniform gap above the ground 
state sector if the spectrum a{H\) of the Hamiltonian H\ satisfies the following conditions: 
The ground state of the Hamiltonian H^. is q-fold (quasi) degenerate in the sense that there 
are q eigenvalues, Eq^i, . . . , Eo^g, in the ground state sector at the bottom of the spectrum 
of H\ such that 

AS := meix{\Eo^^ - Eoy\} ^ as |A,| ^ oo. (2.12) 

Further the distance between the spectrum, {EQi,...,EQq}, of the ground state and the 
rest of the spectrum is larger than a positive constant AE which is independent of the 
volume \As\. Namely there is a spectral gap AE above the ground state sector. 

Let Ax, By be observables with the support X,Y C. Ag, respectively. We say that 
the pair of two observables, Ax and By, is fermionic if they satisfy the anicommutation 
relation, {Ax - By} ~ ior X r\Y — If they satisfy the commutation relation, then we 
call the pair bosonic. 

Define the ground-state expectation as 

(• . := -^Tr (. • .)Po,A, (2.13) 

where Pq.a is the projection onto the ground state sector. For the infinite volume, 

(• • •)o := weak*- hm (• , (2.14) 

|As|Too ' 

where we take a suitable subsequence of finite lattices A going to the infinite volume so 
that the expectation converges to a linear functional for a set of quasilocal observables. 
Although the ground-state expectation thus constructed depends on the STibsequence of 
the lattices and on the observables, our results below hold for any ground-state expectation 
thus constructed. Further, we denote by 

uj(- ■ •) := weak*- lim ($a, (• • O^a) (2.15) 

|As|Too 

the ground-state expectation in the infinite volume for a normalized vector $a in the 
sector of the ground state for finite lattice A. 

Theorem 2.6 (Clustering of fermionic correlations): Let Ax, By be fermionic ob- 
servables with a compact support. Assume that there exists a uniform spectral gap AE > 
above the ground state sector in the spectrum of the Hamiltonian H\ in the sense of Def- 
inition 2.5. Let cu be a ground-state expectation (2.15) in the infinite volume limit. Then 
the following bound is valid: 

u{AxBy) - ^ [u;{AxPoBy)-Uj{ByPoAx)] 



< 



Const X dist(X, Y)] for power-law decaying hx; ,^ -^gx 

[ exp[— /i dist(X, F)], for exponentially decaying hx, 
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where Pq is the projection onto the sector of the infinite-volume ground state, ^ and 

V — ^ , ^ ^ and a — „^ , ^ ^ . (2.17) 

Here and are, respectively, an increasing function of rj and fi, and give an upper 
hound of the group velocity of the information propagation. 

Remcirk: Clearly there exists a maximum /imax such that the bound (2.8) holds for 
any ^ < fi„iax- Combining this observation with (2.17), there exists a maximum fl = 
maXju<^^^^{/x/ (1 + 2Vfj,/AE)} which gives the optimal decay bound. When the interaction 
hx is of finite range, one can take any large But the upper bound V/j, of the group 
velocity exponentially increases as n increases because depends on Aq of (2.8). In 
consequence, a finite jl gives the optimal bound. 

Formally applying the identity, {AxPoBy)o = (-ByPo^x)o, for the bound (2.16), we 
have the following decay bound for the correlation:^ 

Corollciry 2.7 Let Ax, By be fermionic observables with a compact support. Assume that 

there exists a uniform spectral gap AE > above the ground state sector in the spectrum 
of the Hamiltonian H\ in the sense of Definition 2. 5. Then the following bound is valid: 

\{AxBy) I < Const X I ^ dist(X, Y)]'^, for power-law decaying hx; 18) 
' I exp[— /i dist(X, y)], for exponentially decaying hx , 

in the infinite volume limit, where fj, fi are as defined above. 

Theorem 2.8 (Clustering of bosonic correlations): Let Ax, By be bosonic observ- 
ables with a compact support. Assume that there exists a uniform spectral gap AE > 
above the ground state sector in the spectrum of the Hamiltonian H\ in the sense of Def- 
inition 2.5. Let cu be a ground-state expectation (2.15) in the infinite volume limit. Then 
the following bound is valid: 

u{AxBy) - ^ [uj{AxPoBy) + Uj{ByPoAx)] 

< Const X I ^ dist(X, 1")]"^, for power-law decaying hx; ,^ 
~ ' [ exp[— /i dist(X, y)], for exponentially decaying hx , 

where fj, fl are as defined above. 

Remcirk: Theorem 2.8 is a clustering bound for the connected correlation functions. 

We now make some additional definitions that will enable us, in certain cases, to prove 
the decay of [u!{AxPoBy) -\- cu^ByPQAx)] /2 so that Theorem 2.8 can be replaced with a 
stronger bound below. Theorem 2.10. 

^uj{- ■ ■ Pq- ■ ■) is also defined as a bilinear functional for a set of quasilocal observables in the weak* 
limit. 

^See Section 3 for details. 
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Definition 2.9 (Self-similcirity): Write m — with the degeneracy q of the ground 
state sector. We say that the system has self-similarity if the following conditions are 
satisfied: For any observable A of compact support and any given large L > 0, there 
exist transformations, Ri, R2, . . . , Rm, and observables, B'^^\ B^'^\ . . . , B^"^\ such that the 
Hamiltonian H\ is invariant under the transformations, i.e., Rj{H\) = H\ for any lattice 
A with sufficiently large \As\, and that the observables satisfy the following conditions: 

B^^^ = Rj{A) and (b^^))^ = Rj{A'^) for j = 1, 2, . . . , m, (2.20) 

dist(supp A, supp B^^^) >L /or j = 1, 2, . . . , m, (2.21) 

and 

dist(supp S^-'^supp S(^)) > L for j ^ k. (2.22) 

In Section 4, we will discuss other conditions similar to this self-similarity condition. 

Theorem 2.10 Assume that the degeneracy q of the ground state sector of the Hamil- 
tonian H\ is finite in the infinite volume limit, and that there exists a uniform spectral 
gap AE > above the ground state sector in the spectrum of the Hamiltonian H\ in the 
sense of Definition 2.5. Further assume that the system has self- similarity in the sense 
of Definition 2.9, and that there exists a subset Al of bosonic observables with a compact 
support such that Rj{Ai,) G Af, — (Al)^ for j — 1,2, ... ,m, and that {A'xBy)q as 
dist(X, y) — s> 00 for any pair of bosonic observables, A'^.By G Al. Let u be a ground- 
state expectation (2.15) in the infinite volume limit, and let Ax, By he a pair of bosonic 
observables satisfying Ax € Al. Then the following bound is valid: 

\uj{A B )\ < Const X ( ^ dist(X, Y)]'^, for power-law decaying hx; 23) 
' [ exp[— /i dist(X, F)], for exponentially decaying hx, 

where fj, ft are as defined above. 

Remark: 1. If the finite system is translationally invariant in one of the spatial directions 
with a periodic boundary condition, then the self-similarity condition of Definition 2.9 is 
automatically satisfied by taking the translation as the transformation Rj. Thus we do 
not need an additional assumption for such systems. 

2. Theorem 2.10 can be extended to a system having infinite degeneracy of the ground 
state sector in the infinite volume limit if the degeneracy for finite volume is sufficiently 
small compared to the volume of the system. See Theorem 4.1 in Section 4 for details. 

In order to apply this theorem, we need to be able to show that (AxBy)^ — > as 
dist{X,Y) — >• 00 in the infinite volume for any pair of bosonic observables. Ax, By G 
However, this was proven [2] for quantum spin or fermion systems with a global U(l) 
symmetry on a class of lattices with (fractal) dimension D < 2 as defined in (2.25) below, 
so long as the observables behave as a vector under the U(l) rotation. 

We first define the dimension for these lattices. The "sphere", Sr{x), centered at a; £ 
with the radius r is defined as 



Sr{x) {y e A5|dist(y,x) = r}. 



(2.24) 
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Assume that there exists a "(fractal) dimension" D > 1 of the lattice A such that the 
number |5'j.(x)| of the sites in the sphere satisfies 

sup |5^(a;)| < Cor^~^ (2.25) 

with some positive constant Cq. This class of the lattices is the same as in [8]. 

Consider spin or fermion systems with a global U(l) symmetry on the lattice A with 
(fractal) dimension 1 < D < 2, and require the existence of a uniform gap above the ground 
state sector of the Hamiltonian H\ in the sense of Definition 2.5. Although the method 
of [2] can be applied to a wide class of such systems, we consider only two important 
examples, the Heisenbcrg and the Hubbard models. We take the set to be the bosonic 
obscrvables which behave as a vector under the U(l) rotation. In the rest of this section we 
use the results of [2] to show as in (2.29,2.31) that the correlation function for this class of 
observables in these models does decay to zero as dist(X, Y) — > oo. The bounds (2.29,2.31) 
provide only a slow bound on the decay. However, this slow bound on the decay suffices, 
in conjunction with the self-similarity condition of Definition 2.9 to apply Theorem 2.10. 
Thus, under the self- similarity assumption as well as the gap assumption, all the upper 
bounds below (2.29,2.31) are replaced with exponentially decaying bounds by Theorem 2.10. 
In particular, a system with a translational invariance automatically satisfies the self- 
similarity condition as mentioned above. Therefore the corresponding correlations show 
exponential decay for translationally invariant systems on one- dimensional regular lattices. 



XXZ Heisenberg model: The Hamiltonian Ha is given by 

HA^Hr + VA{{Si'^) (2.26) 

with 

{x,y}eAi, 

where = {Si^\ , S^^^) is the spin operator at the site a: e A^ with the spin S = 
1/2, 1,3/2, . . ., and J,^J are real coupling constants; Va{{SI;^^}) is a real function of the 
^-components, {S^^}xeAs^ of the spins. For simplicity, we take 

{x,y}£At 

with real coupling constants J^y. Assume that there are positive constants, and J^ax) 
which satisfy \J^^\ < and \J^^y\ < J^ax bond {x,y} G A;,. 

Consider the transverse spin-spin correlation, (^S^Sy'^^, where S"^ := S*^^-* ±iSj^\ 

Theorem 2.11 Assume that the fractal dimension D of (2.25) satisfies 1 < D < 2, and 
that there exists a uniform spectral gap AE > above the ground state sector in the 
spectrum of the Hamiltonian H\ of (2.26) in the sense of Definition 2.5. Then there exists 
a positive constant 7 such that the transverse spin-spin correlation satisfies the bound, 

< Const, exp f-7{dist(a;, y)Y~'^''^\ , (2.29) 



X y / Q 

in the thermodynamic limit \A.s\ 00. 
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The proof is given in [2], and we remark that the result can be extended to more 



complicated correlations such as the multispin correlation, ^5"+ • • • S^.Sy_^ ■ ■ ■ Sy^J^. If the 
system satisfies the self-similarity condition of Definition 2.9, then the upper bound, (2.29), 
can be replaced with a stronger exponentially decaying one by Theorem 2.10. 

Hubbard model [9, 10]: The Hamiltonian on the lattice A is given by 

= - E E {ixA,aCy,a + tUy^aCx,a) + V {{fl,,^}) + ^ " S„ (2.30) 

where cj, ^, c^^a ai'e, respectively, the electron creation and annihilation operators with the 
z component of the spin ^ =]•[• n^.a = cl:a'^x,a is the corresponding number operator, 
and Sx = {S^^\ , S^^^) are the spin operator given by 5^") = J2a,/3=-\,i4:,a'^a^i3'^x,i3 with 
the Pauli spin matrix (cr^"^) for a = 1, 2, 3; tij e C are the hopping amplitude, V{{nx,a}) 
is a real function of the number operators, and Bj; = (B^\ B^^\ B^^^) e are local 
magnetic fields. Assume that the interaction V{{nx^a}) is of finite range in the sense of 
the graph theoretic distance. 

Theorem 2.12 Assume that the fractal dimension D of (2.25) satisfies 1 < D < 2, 
and that there exists a uniform spectral gap AE > above the ground state sector in 
the spectrum of the Hamiltonian H\ of (2.30) in the sense of Definition 2.5. Then the 
following bound is valid: 



< Const, exp -7{dist(x, i/)}^"^/^ (2.31 



with some constant 7 in the thermodynamic limit \As\ 00. // the local magnetic field 
has the form B^. = (0, 0, B^), then we further have 

S^S;)^ < Const, exp [-y{dist(x, y)}^"^/'] (2.32) 



with some constant 7'. 



The proof is given in [2]. Clearly the Hamiltonian H\ of (2.30) commutes with the 
total number operator A/a = J2x€A-s S/i=T,J. ^ finite volume |As| < 00. We denote 

by H\^N the restriction of i^A onto the eigenspace of A/a with the eigenvalue A^. Let Po.A,7V 
be the projection onto the ground state sector of i?A,Ar, and we denote the ground-state 
expectation by 

(• • ■)o,u = weak*- lim — Tr (• • ■)Po,a,n, (2.33) 

where is the degeneracy of the ground state, and u is the limit of the filling factor 
A^/|As| of the electrons. Since the operators 5"^ do not connect the sectors with the 
different eigenvalues N, we have 

Theorem 2.13 Assume that the fractal dimension D of (2.25) satisfies 1 < D < 2, 
and that there exists a uniform spectral gap AE > above the ground state sector in the 
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spectrum of the Hamiltonian H\ n in- the sense of Definition 2. 5. Then the following hound 
is valid for the filling factor v of the electrons: 



(2.34) 



< Const, exp -7'{dist(a;,y)}^"-^/^ 

with some constant 7' in the infinite volume limit. 

The proof is given in [2]. If the system satisfies the self- similarity condition of Defini- 
tion 2.9, then these three upper bounds, (2.31), (2.32) and (2.34), can be replaced with a 
stronger exponentially decaying one by Theorem 2.10. 



3 Clustering of correlations 

In order to prove the power-law and the exponential clustering, Theorems 2.6 and 2.8, 
we follow the method [3]. The key tools of the proof are Lemma 3.1 below and the Lieb- 
Robinson bound [6, 11] for the group velocity of the information propagation. The sketch 
of the proof is that the static correlation function can be derived from the time-dependent 
correlation function by the lemma, and the large- distance behavior of the time-dependent 
correlation function is estimated by the Lieb-Robinson bound. As a byproduct, we obtain 
the decay bound (2.18) for fermionic observables. 

Consider first the case of the bosonic observables. Let Ax, By be bosonic observables 
with compact supports X,Y G A^, respectively, and let Ax{t) = e'^'^^'^AxC^'*^''^, where 
t e R and H\ is the Hamiltonian for finite volume. Let $ be a normalized vector in the 
ground state sector. The ground state expectation of the commutator is written as 

{^,[Ax{t),By]^) = (<l>,Ax(t)(l-Po,A)5y<l>)-($,fiy(l-Po,A)Ax(t)$) 

-I- Ax {t)Po,ABY^) By Po^AAxm) . (3.1) 

In terms of the ground state vectors $o,i^, = 1, 2, . . . , with the energy eigenvalues, Eq^^, 
and the excited state vectors with En, n = 1,2, . . ., one has 

($, Ax{t){l - Po,A)5y$) = E E (^0,., Ax^n) ($n, By^oy) e'^^^^--^^-) , (3.2) 
($, By{l - Po,A)Axm) = E E (^0,., 5y$n) ($n, Ax<^oy) e**(^"-^o-'), (3.3) 

Ax{t)Po,ABy^) = E E (^0,., Ax^o,^) ($o,M, By^oy) e-^*(^°--^°-) (3.4) 
vy M 

and 

ByPo,AAxm) = E E (*o,., By^o,i.) ($0,^, Ax^oy) e^*(^°--^o,.')^ (3.5) 
i^y fj- 

where we have written ^ 

$ = Ea.*o,.- (3.6) 

In order to get the bound for {^,Ax{t — 0)5y$), we want to extract only the "neg- 
ative frequency part" (3.2) from the time-dependent correlation functions (3.1). For this 
purpose, we use the following lemma [3]: 
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Lemma 3.1 Let E eTi, and a > 0. Then 

1 [W /-o 



lim lim — 



T ^—iEt^—at'' 



T t + ie 



dt 



f l + C>(exp[-AEV(4a)]) forE>AE; 
\ C(exp[-AEV(4a)]) for E < -AE. 



Proof: Write 

m 

Using the Fourier transformation, 



2^ 



T ^—lEt^—at^ 

-T t + ie 



dt. 



r exp[-(u; + E)y{Aa)]e'^'du;, 



we decompose the integral I{E) into three parts as 

I{E)^l4E) + IoiE) + I+{E), 

where 



I4E) 



i I [W fT ^ 1 
dt- 



27r 27r V q; 



t + ie J-oo 



-Auj 



du;ex.p[-{u; + Ef/{Aa)]e 



iu!t 



and 



' 27r27rVQ; 



+ ^ ^ 27r27rVQ; 



dt 

T t + 



-T 



dt 

■T t + 



— / da;exp[-(u; + E)7(4a)]e 

h le J-Au! 

-I 

h le 



dujexp[-{<jj + EY/{Aa)y\ 



where we choose Ao; = hT ^1"^ with some positive constant h. 
First let us estimate Io{E). Note that 



ze 



t + ie t^ + e^ t^ + e2 



Using this identity, one has 

i 1 fW /-^^ 



7 1 pjF rAu) pT 



t sin cut ie cos Lvt 



t2 + e2 



t2 + e2 



(3.7) 
(3.8) 
(3.9) 
(3.10) 

(3.11) 

(3.12) 

(3.13) 



(3.14) 



(3.15) 



where we have interchanged the order of the double integral by relying on \t\ < T < oo. 
Since the integral about t can be bounded by some constant, one obtains 



|7o(^)| < Const. X a-^/^Au; < Const, x a-^/^T"^/^. 

Therefore the corresponding contribution is vanishing in the limit T | cxo. 
Note that 



l!L = / C>(cu-'T-i) for > 0; 

2n J-T ""t + ie~ \ e''' + 0{uj-^T-^) for cu <Q. 



— / dt- 



(3.16) 



(3.17) 
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Using this, the function I+{E) of (3.13) can be evaluated as 

|/+(^)| < Const. X r-^/2_ 



(3.18) 



This is also vanishing in the limit. 

Thus it is enough to consider only the integral I-{E). In the same way as the above, 
one has 

I_{E) = — J- / cia;exp[-(u; + Ef/{Aa)]e'^ + 0{T-^'''). (3.19) 

ZTT V Q; J-oo 

Since e^'^ < 1 for a; < 0, one has 

limlim7_(£;) = ^y^ r duex^[-{u + Ef / {Aa)]. (3.20) 

Ttoo ej,0 M \ a J-oo 



Note that, for E < -AE, 

f dueM-i^ + ^)V(4«)] < ^exp[-AEV(4a)], 

ZTT V Q; J-oo 2 

and, for E > AE, 

-!-,/-/ cia;exp[-(u; + £;)V(4a)] = 77-^/-/ (iu;exp[-(u; + £;)V(4q;)] 
27r V q; J-oc 2tt \ a J-oo 

1 ^tT 

— J- duexp[-{u; + Ef/{4a)] 



(3.21) 



27r V a 

= l + C(exp[-AE7(4a)]). 



(3.22) 



Clearly these imply (3.7). 



From Lemma 3.1 and the expression (3.1) of the correlation function with (3.2) and 
(3.3), one has 



dt- 



1 



($,[Ax(i),Sy]$) e 



limlim— , 

TToo eio 2ti J-t t + ie 

= ($, Ax{l - Po,a)By^) + 0(cxp[-A^V(4a)]) 

+ lim lim ^ r dt-^ [($, Ax{t)Po,ABY^) - 5yPo,A^x(i)$)] e-"*'(3.23) 

Ttoo eiO ZTT J-T t + l€ 

for finite volume. 

In the following, we treat only the power-law decaying interaction hx because one can 
treat the exponentially decaying interactions in the same way. See also refs. [3, 4] in which 
the exponential clustering of the correlations is proved for finite-range interactions under 

the gap assumption along the same line as below. 

In order to estimate the left-hand side, we recall the Lieb- Robinson estimate (A.l) in 
Appendix A, 



< Const. X 



3^1*1 _ I 



(1 + r)'? \t\ 



(3.24) 
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for r > 0, where we have written r = dist(X, Y). Using this estimate, the integral can be 
evaluated as 



dt 



< 



t-\- ie 
[^,[Ax{t),By]^) 



dt 



'\t\<ce 
< Const. X 



t + ie 



+ 



1 



+ 



Const. 



exp[— ac 



\t\>c£ 



^^ {<^,[Axit),BY]<^>) ^_^,. 
t + ie 



where c is a positive, small parameter, and i — log(l + r), and we have used 

e"l*l - 1 



(3.25) 



\t\<c£ \t\ 



-dt < 2e 



(3.26) 



In order to estimate the integral in the right-hand side of (3.23), we consider the matrix 
element ($o,i/, ^x(^)-fo,A-By$o,i^') because the other matrix elements in the ground state 
can be treated in the same way. Using Lemma 3.1, one has 



hm hm ^ r dt-^ ($0,., Ax{t)Po,ABY^oy) e""*' 
rtoo eio 27r J-T t + ie 

r dcuexp[-(cu + A£^,,)V(4a)], 

J —oo 



1 1 /tT 



27r V a 



(3.27) 



where AS^^i, = Eq^^ — Eq^^. Using the assumption (2.12) and the dominated convergence 
theorem, we have that, for any given £ > 0, there exists a sufficiently large volume of the 
lattice A, such that 



lim lim 



T 



1 



dt—- {%^^,Ax{t)Po,ABY%y) - - {^o,^,AxPo,aBy^o,^ 



< e. (3.28) 



T^ooelO 27r J-T t + ie' ' ' 2 

Combining this observation, (3.23) and (3.25), and choosing a = AE/{2ci), one obtains 



u;{AxBy) - \ \uj{AxP^By) + uj{ByP^Ax)\ 



< Const. X 



[1 + r)^- 



+ Const. X exp 



cAE 



(3.29) 



in the infinite volume limit, where the ground-state expectation uu is given by (2.15). 
Choosing c = ri/{v + AE/2), we have 



uj{AxBy) - - [uj{AxPoBy) + uj{ByPoAx)] 



with fj — r]/{l + 2v/ AE). In the same way, we have 



< 



Const. 



[1 + dist(X, Y)Yi' 



(3.30) 



uj{AxBy) - - [uj{AxPoBy)+uj{ByPoAx) 



< Const. X exp[-/i dist(X, ¥)] (3.31) 
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for the exponentially decaying interaction hx, where jl — ///(I + 2v/AE). This proves 
Theorem 2.8. The corresponding bound for finite-range interactions was already obtained 
in [4]. Using the definition (2.13) of the expectation (• • •)o,a and the identity, 

(AxPoBy),^^ = {ByPoAx)o,a , (3-32) 
for the integral in the right-hand side of (3.23), we obtain 

< Const X / ^ dist(X, Y)]-^, for power-law decaying hx] /g 
~ ' I exp[— /i dist(X, y)], for exponentially decaying /ix 

for any finite lattice As D X,Y in the same way as in the above. 

Next consider the case that the pair, Ax, By, is fermionic. Note that 

{%,,,{Ax(t).By}<^0..') 

= ($0,., Axit)il - Po.A)BY^o.y) + ($0,., By{1 - Po,A)Ax{t)^oy) 

+ ($0,., Ax{t)Po,ABY^oy) + (*o,., BYPo,AAx{t)^oy) ■ (3.34) 

Since the difference between bosonic and fermionic observables is in the signs of some 
terms, one has 

uj{AxBy) - ^ [uj{AxPoBy) - uj{ByPoAx)] 

< Const X I ^ dist(X, Y)]'^, for power-law decaying hx] /g 
~ ' [ exp[— /i dist(X, y)], for exponentially decaying /ix- 

In particular, thanks to the identity (3.32), we obtain 

\(A B ) \ < Const X I ■^''l"''' for power-law decaying /ix; 

° ~ ' [ exp[— /i dist(X, y)], for exponentially decaying /ix- 

(3.36) 

This is nothing but the desired bound. We stress that, for infinite degeneracy of the 
infinite-volume ground state, this upper bound is also justified in the same argument with 
the dominated convergence theorem. 



4 Vanishing of the matrix elements in the ground 
state 

The aim of this section is to prove the bound (2.23) for the correlation and discuss an 
extension of Theorem 2.10 to a system having infinite degeneracy of infinite- volume ground 
state. The latter result is summarized as Theorem 4.1 below. We will give only the proof 
of Theorem 4.1 because Theorem 2.10 is proved in the same way. By the clustering bounds 
(3.30) and (3.31), it is sufficient to show that all the matrix elements, {^oyAx^o,u), in 
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the sector of the ground state are vanishing. The key idea of the proof is to estimate the 
absolute values of the matrix elements by using the self-similarity condition and the decay 
bound (4.2) below of the correlations at a sufficiently large distance. 

Wc denote by the degeneracy of the sector of the ground state for the finite lattice 
A, and wc allow q\ ^ oo as \As\ | oo. We write m = q\ for short. To begin with, we 
write the bound (3.33) as 

(^xSy)o,A - (^xPo,A5y)o,A| < Go(dist(X, Y)), (4.1) 

where we have written the upper bound of the right-hand side by the function Gq of the 
distance. We assume that the following bound holds: 

{AxBy)o^^\ < G'i(dist(X, Y)) (4.2) 

with an upper bound Gi which is vanishing at the infinite distance. Further we define Ga 
as 

Ga{Ax,By) :=max{G'o(dist(X,r)),G'i(dist(X,r))}. (4.3) 

Theorem 4.1 Let cu be a ground-state expectation (2.15) in the infinite volume limit, and 
let Ax, By he a pair of bosonic observable with compact supports X, Y. Assume that there 

exists a uniform spectral gap AE > above the ground state sector in the spectrum of the 
Hamiltonian in the sense of Definition 2.5. Suppose that, for any given e > 0, there 
exists Mq > such that, for any large lattice A satisfying \As\ > Mq, there exists a set 
of observables, B^\ j — 1, 2, . . . , m, and a set of transformations, Rj, j = 1, 2, . . . , m, 
satisfying the following conditions: Any pair of the observables. Ax , B^^\ . . . , B^"^\ is 
bosonic, 

B^^^^Rj{A), {B^^'^y ^ Rj{A^x) and Rj{Hx) = H^, (4.4) 

and 

where we have written B^^^ — Ax. Then we have the bound, 

\uj(AxBy)\ < Const x | ^ dist(X, Y)]-^, for power-law decaying hx; g^ 
' I cxp[— /i dist(X, F)], for exponentially decaying hx , 

in the infinite volume limit. 

Proof: From the bound (3.30) or (3.31) and the Schwarz inequality, 

\u;{AxPoBy)\' < u;{AxPoA^xMbIPoBy), (4.7) 
it is sufficient to show uj{AxPoAx) — uj{AxPoAx) — 0. Further, we have 

$, APo,aA^^) < qA {APo,aA^)^^^ = qA {A^Po,aA)^^^ (4.8) 
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for any ground state vector $ with norm one and any observable A on the finite lattice A. 
Therefore we estimate (^x-Po.a^x)^^- 

Note that, from the clustering bound (4.1), (4.2) and (4.3), we have 



{AxPoaBy). 



< 



(AxBy) 



0,A 



0,A 



We define 



< 2G^{Ax,By). 



(4.9) 



(4.10) 



for J = 0, 1, . . . , m and for the finite lattice A, where we have written i — {i/', u) with 
i — l,2,...,m for short. Since {Bi\ ■ ■ ■ , B^'') is an m-dimensional vector, there 
exist complex numbers, Cj,j — 0, 1, . . . , m, such that, at least, one of Cj is nonvanishing 
and that 

m 
j=0 

Let i be the index which satisfies jC^I = max{|Co|, |Ci|, . . . , |Cm|}- Clearly, we have 



(4.11) 



(4.12) 



Therefore 



-t m „ 

{B^'^yPoAB^'^ =-E^P 



f^e qA i=i 

{B^'^yPoAB^'^ 



< mmax 



0,A 



< 



2glmax{GA((SW)t,sO))}, (4.13) 



where we have used the inequality (4.9) for getting the last bound. When £ = 0, we obtain 

?a(AxPo,a^^)„ <2e (4.14) 



0,A 



from fi*^*^^ = Ax assumption (4.5). When £ 0, we reach the same conclusion by 

using the relation. 



^xPoaAx)^^ = {RM^x)PoaRM))^^ = nB('^yPo,AB 



0,A 



(4.15) 



which is derived from the assumption (4.4). 



Remark: 1. The advantage of Theorem 4.1 is that it is easier to find and Rj because 
of the finiteness of the lattice. Actually one can construct B^^\ Rj and A satisfying 
the requirement by connecting m copies of a small, finite lattice to each other at their 
boundaries. But, if the degeneracy q\ exceeds y''|As|, we cannot find the observables, B^^\ 
and the transformations, Rj. Therefore our argument does not work in such cases. 
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2. Under the weaker assumption, 



q\ max 

i,je{0,l,...,m}: 



{ 



G^((S«)^E(^■))}<6, 



(4.16) 



than (4.5), we can obtain the bound. 



\{AxBy)q\ < Const. X 



( 



[1 + dist(X, Y)] for power-law decaying hx', 
exp[— /i dist(X, Y)], for exponentially decaying hx, 



(4.17) 



in the infinite volume limit. 

3. Consider the situation of the above Remark 2 or the case with a finite degeneracy of 
the infinite- volume ground state. Then, instead of introducing the transformations Rj, we 
can directly require 



in the infinite volume limit, and at infinite distance between the observables Ax and 



A Lieb-Robinson bound for group velocity 

Quite recently, Nachtergaele and Sims [6] have extended the Lieb-Robinson bound [11] 
to a wide class of models with long-range, exponentially decaying interactions. In this 
appendix, we further extend the bound to the power-law decaying interactions. We also 
tighten the bound on the exponentially decaying case. (See Assumption 2.4 compared to 
that in [6]). However, in our proof, the time t must be real. 

In the following, we treat only the case with bosonic observables and with the power- 
law decaying interaction hx because the other cases including the previous results can be 
treated in the same way. 

Theorem A.l Let Ax, By be a pair of bosonic observables with the compact support, 
X,Y, respectively. Assume that the system satisfies the conditions in Assumption 2.1 or 



where the positive constants, C and v, depend only on the interaction of the Hamiltonian 
and the metric of the lattice. 

Remcirk: The same bound for fermionic observables is obtained by replacing the com- 
mutator with the anticommutator in the left-hand side. 




(4.18) 



2.2. Then 



\\[Ax{t),BY]\\<C\\Ax\\\\BY\\\X\\Y 



[1 + dist(X, Y)]n 



for dist(X,y) > 0, (A.l) 



For exponentially decaying interaction hx, the following bound is valid: 
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Theorem A. 2 Let Ax, By he a pair of bosonic observables with the compact support, 
X,Y, respectively. Assume that the system satisfies the conditions in Assumption 2.3 or 
2.4. Then 

\\[Ax{t),BY]\\ < C\\Ax\\\\BY\\\X\\Y\exp[-ijdist{X,Y)] [e'''*! - l] for dist(X,F) > 0, 

(A.2) 

where the positive constants, C and v, depend only on the interaction of the Hamiltonian 
and the metric of the lattice. 

Reiricirk: For the proof under Assumption 2.3, we rely on the inequahties, (2.9) and 
(2.10). Assumption 2.4 is milder than that in ref. [6] as remarked in Section 2. 

We assume that the volume | A^l of the lattice A is finite. If it is necessary to consider the 
infinite volume limit, we take the limit after deriving the desired Lieb- Robinson bounds 
which hold uniformly in the size of the lattice. Let A, B be observables supported by 

compact sets, X,Y C A^, respectively. The time evolution of A is given by A{t) = 
^itHAj^^-ttHA _ YiYst, let us derive the inequality (A. 12) below for the commutator [A(t), B]. 
We assume t > because the negative t can be treated in the same way. Let e = t/N 
with a large positive integer N, and let 

t^^Ln for n = 0, 1,...,A^. (A.3) 



Then we have 



mt),B]\\ - ii[A(o,. Bill = s . X mtM\\-\mK)M\ ,^ 



i=0 



In order to obtain the bound (A. 12) below, we want to estimate the summand in the 
right-hand side. To begin with, we note that the identity, ||C/*0[/|| = ||0||, holds for any 
observable O and for any unitary operator U . Using this fact, we have 

mtn^ilB]\\ - mt^lB]\\ = \\[A{e),B{-tn)]\\ - \\[A,B{-t^)]\\ 

< \\[A + te[Hj,, A] , B{-t^)] II - II [A, i?(-U)] || + 0(6^) 
= \\[A + te[Ix,A],Bi-tr.)]\\ - \\[A,Bi-tr.)]\\ + 0{e') 



with 



where we have used 



(A.5) 



E hz, (A.6) 

z-.znxj^^ 



A{e)^A + ie[HA,A] + 0{e'') (A.7) 
and the triangle inequality. Further, by using 

A + ie[Ix, A] = e^^-^^Ae"^^^^ + ^(e^), (A.8) 
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we have 



\\[A + ie[Ix,AiB{-t^)]\\ < \[e'^'-Ae-'''-,B{-tn)]\\ + 0{e') 

[A,e-'''^B{-tn)e'''^] +0(6") 



< \\[A, B{-ti) - ie[Ix, Bi-t„)]] II + Oie') 

< ||[ASH„)]||+e||[A[/x,SH„)]]||+0(e^). (A.9) 



Substituting this into the right-hand side in the last hne of (A. 5), we obtain 

||KWi),5]||-||[^(^n),S]|| < e\\[A[Ix,B{-t^m + 0{e') 

< 26||A||||[/^(t„),S]|| + 0(e^). 



(A.IO) 



Further, substituting this into the right-hand side of (A. 4) and using (A. 6), we have 



N-l 



||[A(t),i?]||-||[A(0),i?]|| < 2p||Eex||[/x(tn),i?]|| + 0(6) 



n=0 



N-l 



< 2\\A\\ J2 J2^^\\[hz{tn),B]\\ + 0{e). (A.ll) 

Z:ZnX^0 n=0 

Since hz{t) is the continuous function of the time t for a finite volume, the sum in the 
right-hand side converges to the integral in the hmit e J, (A^ t oo) for any fixed finite 
lattice A. In consequence, we obtain 



\\[A{t),B]\\ - \\[A{0),B]\\ < 2\\A\\ rds\\[hz{s),B]\\. 



(A.12) 



We define 



CB{X,t) := sup 



\mt),B]\\ 



(A.13) 



where Ax is the set of observables supported by the compact set X. Then we have^ 

r\t\ 

\\hz\\ I 



CB(X,i) <Cb(X,0) + 2 W^zW dsCB{Z,s) (A.14) 



from the above bound (A.12). 

We recall that the observables, A and B, are, respectively, supported by the compact 
sets, X, F C A^. Assume dist(X, Y) > 0. Then we have Cb{X, 0) = from the definition 
of Cb{X, t), and note that 



Cb(Z,0) < 



2||5||, forZnr^0; 
0, otherwise. 



(A.15) 



Since the loeal interaction hz with Z C X does not change the support X of ^ in the time evolution, 
we can expect that the sum in the right-hand side of (A.14) can be restricted to the set Z satisfying 
Z n X 7^ and Z \ X 7^ 0. However, this restriction does not affect the resulting Lieb- Robinson bound. 
Therefore we omit the discussion. 
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Using these facts and the above bound (A. 14) iteratively, we obtain 

\\hzA\ / 

Zi:ZinXj^9 
Zv.Zinx^fD 

.\t\ H^il 



CB{X,t) < 2 Yl \\hzA\ rdsiCBiZi,si) 



r\t\ r\si\ 

+ 2^ E ll^^lll E 11^2211/ dsj dS2CBiZ2,S2) 
Zi:ZinX^0 Z2:Z2nZi^(D ''^ 

< 2||B||(2|t|) E ll^^ill 

+ 2||i^||^ E II^^J E II^^J 

Zr-ZiDX^H) Z2;Z2nZi7^0,Z2nY7^0 

+ 2||s||^ E \\hzA\ E 11^^.11 E 11^^311 + 

Zv-ZinXj^$ Z2:Z2nZi^9 Zz-.Z^^Z^+'h^ZzV^Y^^ 



(A.16) 



Proof of Theorem A.l under Assumption 2.1: The first sum in the power series (A.16) is 
estimated as 

E \\hzA\ < \\hzA 

Zx-.Zx nX^0,Zi ny 7^0 x&y& Zi 3x,y 

- [1 + dtf('x^y)]'^ ^^-^^^ 

from the assumption (2.2). The second, double sum is estimated as 

E ll^^ill E 11^22 II 

Zx:Zx<r\X+% Z2:Z2r\Zi^%,Z2nY^% 

< E \\hzA E 11^^211 

xeX J/gy Zl2 6As Z\3x,Z\2 Z23zi2,y 

" kkhz^^K [1 + ^12)]" [1 + dist(^i2, y)]"^ 

- [i + dist(x,y)]'^' ^^-^^^ 

where we have used the assumptions (2.2) and (2.3). Similarly, the third, triple sum can 
be estimated as 

E ll^^ill E 11^^2 11 E ll^^sll 

Zi:Zinx^0 Z2:Z2nZi^0 Z3:Z3nZ27^0,Z3ny7^0 

< EE E E E II^^JI E 11^^211 E ll^^all 

xeX yeY Zl2&h.s Z23eK Zi3x,Zl2 Z2BZ12,Z23 Z33Z23,y 
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^tx z,%s [1 + dist(x, Zi2)]^ [1 + dist(^i2, Z2z)]^ [l + dist(^23, V)]^ 



< V V V 

~ iiyi^.i^A. [1 + dist(x,zi2)]'' [1 + dist(zi2,y)]'' 



- [i + dist(x,y)]'?- ^'^■'''^ 

Prom these observations, we have 

^^^^''^ - [1 + dist(X, Y)]^ r'^l^° + -^\Po + ^AoPo + • • - j 

2poi5|||X||F| 



[1 + dist(X,F)]'' 
Consequently, we obtain 



{exp[2AoPo|i|] -1}. (A.20) 



I'f^^^^'^^l' - ^Ii?diil!pr,y)ir' WPAoPo|t|] - 1} (A.21) 

from (A. 13). ■ 

Proof of Theorem A.l under Assumption 2.2: The first sum in the power series (A. 16) is 
estimated as 

E 11^^.11 < EE E \\hzA 

Zv.ZinX^0,ZinYj^lll x€X y€Y ZiBx,y 

< EE E ||/izJ|[l + dist(x,y)]-''[l + diam(ZOr 

xeX yeY Zi3x,y 

< [l + dist(X,F)]-''|X||F|so, (A.22) 

where 

So = sup ^ \\hz\\[l + diam(Z)]''. (A.23) 

^ ZBx 

Clearly, this constant sq is finite from the assumption (2.7). The second, double sum is 
estimated as 

E ll^^ill E 11^^2 11 

Zi:Zinx^0 Z2:Z2nZi^0,Z2ny^0 

< EE E E II^^JI E 11^^.11 

xeX yeY 2:12 eAs Zi3x,zi2 Z23zi2,y 

< EE E [l + dist(a;,zi2)]-''[l + dist(^i2,y)]-^ 

xeX yeY zrzeAs 

X E ||/izJ|[l + diam(Zi)f ^ ||/iz2||[l + diam(Z2)]^ 

Zi3x,Zl2 Z2BZl2,y 

< [1 + dist{X, Y)rY.i: E E ||/^zJ|[l + diam(ZOr 

xeX yeY zi2eAs Zi3x,zi2 
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X E ||/iz.||[l + diam(Z2)r 

Z23zi2,y 

< [l + <hst{X,Y)]-^Y.i: E \\hzM^ + di^MZi)]' 

xGX y£Y ZiBx 

X E ll^22lll^2|[l + diam(Z2)]'' 
Z23y 

< [l + dist{X,Y)]-'^\X\\Y\soSi, (A.24) 

where we have used the assumption (2.7) and the inequahty, 

[1 + dist(a;, z)]-'^[l + dist(z, I/)]-" 
= [1 + dist(a;, z) + dist(2;, y) + dist(a;, z)dist{z, y)]"^ 

< [1 + dist(x, + dist(^,|/)]"'' 

< [l + dist(x,y)]-'', (A.25) 
for any z & Ag. Similarly, the third, triple sum can be estimated as 

E ll^^ill E W^z^w E ll^^sll 

Zi:ZinX^0 Z2:Z2nZi^0 Z3:Z3nZ2it<D,Z3nY^<D 

< EE E E E II^^JI E 11^^.11 E ll^^all 

xex yeY zi2eAs z23eAs Zi3x,zi2 -^23^12,223 Z3Bz23,y 

< [l + dist(X,F)]-^;^^ E E E \\hzM^ + diMZi)r 

xGX yeV zi2eAs 223eAs Zi3x,zi2 

X E \\hz2\\[l + dia.m{Z2)r E ||/iZ3l|[l + diam(Z3)P 

•^2 3^12 ,223 Z33Z23,y 

< [l + dist{X,Y)]-'^\X\\Y\sosj. (A.26) 
Prom these observations, we have 

2sosi'\\B\\\X\\Y\ ^ i2s,\t\r 
^^^^''^ - [l + dist(X,F)]'^^^,^^ 

2sos{'^\\B\\\X\\Y 



[1 + dist(X,F)]'? 

As a result, we obtain 



{exp[2si|i|] - 1}. (A.27) 



\\[m,B]\\ < {exp[2.,|t|] - 1} (A.28) 



[i + dist(x,y)]^ 

from (A. 13). ■ 
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